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Abstract 

We calculate the squeezed back-to-back correlation (BBC) functions of (jxj) and K^K~ for heavy- 
ion collisions at RHIC and LHC energies, using (2-|-l)-dimensional hydrodynamics with fluctuating 
initial conditions. The BBC functions averaged over event-by-event calculations for many events 
for the hydrodynamic sources are smoothed as a function of the particle momentum. For heavy-ion 
collisions of Au-|-Au at = 200 GeV, the BBC functions are larger than those for collisions of 

Pb-|-Pb at ylsfpM = 2.76 TeV. The BBC of fxj) may possibly be observed in peripheral collisions at 
the RHIC and LHC energies. It is large for the smaller sources of Cu-|-Cu collisions at = 200 

GeV. 
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I. INTRODUCTION 


In the late 1990s, it was shown [l, 2| that the mass modihcation of the particles in the hot 
and dense hadronic sonrces can lead to a squeezed back-to-back correlation (BBC) of boson- 
antiboson pairs in high-energy heavy-ion collisions. This BBC caused by the interactions in 
medium is different from the pure quantum statistical correlations between the bosons with 
different isospins j^, which are negligible in high-energy heavy-ion collisions because the 
J^ourier transformation of source-density, |/i(0, m)|, is very small even for the lightest meson 
, 1^. Since it is associated with the source medium, the investigations of the BBC may 
possibly provide another way for people to understand the thermal and dynamical properties 
of the hadronic sources formed in high-energy heavy-ion collisions, in addition to particle 
yields and spectra. 


The BBC function is dehned as 


M 


C{k, -k) = 1 + 


|G,(k, -k)p 


a(k,k)G,(-k, -k)’ 

where Gc(ki,k 2 ) and Gs(ki,k 2 ) are the chaotic and squeezed amplitudes, respectively, 

Gc(ki,k2) = (aici^ka), 

Gs(ki,k2) = v'Wkphl^ (OkiOka), 


( 1 ) 


( 2 ) 


( 3 ) 


where Ok and aj^ are the annihilation and creation operators of the free boson with momen¬ 
tum k and mass m, (• • ■) indicates the ensemble average. For a homogeneous source with 
volume V and temperature T, the BBC function can be written as 

V\ci^slnu + c_k s*_y. (n_k + 1)P 


^(k, -k) = 1 + 


( 4 ) 


V [ni(k)ni(-k)] 

where Ck and Sk are the coefficients of Bogoliubov transformation between the creation 
(annihilation) operators of the quasiparticle in medium with modihed mass m* and the free 
observed particle [l, 2|, is boson distribution, 

1 


rik = 


(Uk = + 


and 


exp(Uk/F) - 1’ 
ni(k) = Ickl^Uk + |s_k|^(n_k + 1). 


( 5 ) 


( 6 ) 
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The BBC function ^(k, —k) will be 1 if there is no mass modihcation. However, for a hnite 
mass modihcation, = rri^ — 'ml, /k ~ (5mV(4k^) as |k| —)■ oo, and the BBC function will 


increase with increasing particle momentum [2(|, (^(k, —k) rsj 1+1/k- 


l + kV(^rn74)2. 

the 


For hydrodynamic sources, with the formula derived W Makhlin and Sinyukov j^, 
chaotic and squeezed amplitudes can be expressed as { 2 , 6-^ 


G'e(ki,k 2 ) = 




3^91,2-r 


TT 


■' Xr ' \ e ' ' L-' 


1'^—k^,—k^l ki,—k; 


/ + 
2 


(7) 


G,(k,.k,) = j 


2-7’J J 

o 1^/ 1^/ Ci^/ 1^/ 


X'/l_i,/ ],/ Cl,/ _!,/ 5 _t,/ 1,/ _!,/ “h 1] [*. 

Ki,K2 K 2 K2,k;^l K 2 JJ 


( 8 ) 


Here d^a^{r) is the four-dimension element of freeze-out hypersurface, d^i 2 = 

(fcf -|- ^ 2 ) 72 , and k' is the local-frame momentum corresponding to k* {i = 1,2). The local 


quantities, and n[., are the coefficients of Bogo^ 


\c' Ir' 
Ki,K 2 


iubov transformation and 


boson distribution associated with the particle pair (see in Ref. |8|). 

In Eq. dH]), the factor is equal to for ki = k, k 2 = —ki = —k. So the 

BBC function Cfk, —k) is sensitive to the temporal distribution of the freeze-out points 

iriQ Ti 

2|, l6|-l9[ . Recent research [8| indicates that the BBC functions for the hydrodynamic sources 
with Gaussian initial-energy distributions exhibit oscillations as a function of the particle 
momentum because of the sharp falls of temporal freeze-out distributions at long times. 
Investigating the BBC behavior and predicting its effect in high-energy heavy-ion collisions 
based on more realistic models is of great interest In this work, we investigate the BBC 
functions for the hydrodynamic sources with event-by-event fluctuating initial conditions 


(FIC). Weuse (2+lH 


imensional hydrodynamics with the HIJING 


d FIC and the equation 


12| . and investigate the 


of state s95p-PCE [llj to describe the source evolution as in Ref. 

BBC functions of (fxj) and K^K~ for heavy-ion collisions of Au-|-Au at -^snn = 200 GeV at 
the Relativistic Heavy Ion Collider (RHIC) and Pb-|-Pb at y/sNN = 2.76 TeV at the Large 
Hadron Collider (LHC), respectively. Our investigation indicates that the BBC functions 
averaged over event-by-event calculations for many events for the hydrodynamic sources 
with the FIC are smoothed and without the oscillations appearing in the BBC functions 
for the hydrodynamic sources with Gaussian initial-energy distributions js]. For heavy-ion 
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collisions at the RHIC and LHC energies, the BBC of 00 may possibly be observed in 
peripheral collisions. The investigations for the BBC functions of 00 for Cu+Cu collisions 
at ^/sNN = 200 GeV indicate that the BBC is large for the smaller sources. 

The rest of this paper is organized as follows. In Sec. II, we give a brief review on the 
relativistic hydrodynamic model used in this work, and discuss the space-time distributions 
of particle freeze-out points for the hydrodynamic sources with the FIC. In Sec. Ill, we 
present the BBC functions of 00 and K^K~ for heavy-ion collisions of Au-|-Au and Pb-|-Pb 
at the RHIC and LHC energies, respectively. The dependences of the BBC functions on 
the azimuthal angle and pseudorapidity of the particles are also investigated in this section. 
In Sec. IV, we investigate the source space-time distributions and the BBC functions of 0 
meson for the smaller sources in Cu-|-Cu collisions at ^snn = 200 GeV. Finally, a summary 
and conclusions of this paper are given in Sec. V. 


II. (2 + l)-DIMENSIONAL HYDRODYNAMIC SOURCES 


Relativistic hydrodynamics has been extensively applied in high-energy heavy-ion colli¬ 
sions. In this work, we use the ideal relativistic hydrodynamics in (2-1-1) dimensions to 
describe the transverse expansion of the particle-emitting sources formed in ultrarelativistic 
heavy-ion collisions, and adopt the Bjorken boost-invariant hypothesis 13|] for the longi¬ 
tudinal evolution of the sources. The hydrodynamic equations of motion for the sources 
with zero net baryon density are from the local conservation of energy-momentum Q,Q. 
Under the assumption of Bjorken longitudinal boost invariance, we need only to solve the 
transverse equations of motion in the z = 0 plane 12|, and the hydrodynamic solutions at 
z 0 iy^ = z/t) can be obtained by the longitudinal boost invariance [16|, Il7 |. 

Assuming the local equilibrium of system is reached at time tq, we construct the initial 


energy density of the hydrodynamic source at z = 0, by using t 
the HIJING is used for generating the initial conditions, as [17 


le AMPT code 18 1 in which 
19| 


e(ro,x,2/;z = 0) = AT ^ 


PZc 


Tq 27rcT| 


2 exp 


[x - x^{TQ)f + [y-ya{rQ)\^ 


(9) 


Here is the transverse momentum of parton a in the fluid element at {x,y), Xairo) and 
2/a('ro) are the transverse coordinates of the parton at tq, (Tq is a transverse width parameter, 
and K is a scale factor which can be adjusted to £t the experimental data of produced 
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hadrons [w|. The initial velocity of the flnid element is then determined by the initial 
energy density and the average transverse momentnm of the partons in the element. 

As in Ref 1^, we solve the hydrodynamic eqnations nnmerically by the HLLE scheme and 


the Sod’s operation splitting method 


[l4, 20 


25| . In the calculations, we use the equation of 


state s95p-PCE ll|, which combines the hadron resonance gas at low temperatures and the 


lattice QCD results at high temperatures, and take the parameters cxo = 0.6 fm, tq = 0.6 
and 0.4 fm/c for the heavy-ion collisions at the RHIC and LHC, respectively. With these 
parameter values the hydrodynamic results of transverse momentum spectra and elliptic 
flow of identical pion and charged hadrons are consistent with the experimental dada at the 
RHIC and LHC 


12 


19|. 


For hydrodynamic sources with a Bjorken cylinder, the four-dimension element of freeze- 
out hypersurface can be written as 


(.(’■) = /„(t, ri, 'll) dT(thxd'ii, 


( 10 ) 


where r, r^, and r] are the proper time, transverse coordinate, and space-time rapidity of the 
element. The function /^(t, r^, rj) is related to the freeze-out mechanism that is considered, 
and v) corresponds to the source distributions of proper time and space in the 

calculations [see Eqs. ([7]) and dS])]. In this work we assume that 0 and K mesons are frozen 
out at hxed temperatures Tf = 140 and 160 MeV, respectively, and use the AZHYDRO 


technique [15 



261, 


27| to calculate the freeze-out hypersurface element. 


We plot in Fig. [T] the space-time distributions of the freeze-out points, rx, ?]), 

of (j) meson in the z = 0 plane for central and peripheral heavy-ion collisions of Au-|-Au at 
^/sNN = 200 GeV at the RHIC and Pb-|-Pb at ^/sNN = 2.76 TeV at the LHC. Here, the event 
number for both the RHIC and LHC collisions is 2000. The regions of impact parameter 
b for the denoted centralities 0-10% and 40-50% are taken to be 0-4.2 and 9.0-10.2 fm 
, respectively. One can see that the width of the distributions increases with decreasing 


collision centrality and increasing collision energy. For the hydrodynamic sources with FIC, 
the freeze-out distributions are more dispersive compared to the hydrodynamic sources with 
the Gaussian initial conditions jsj. In Fig. |2l we plot the space-time distributions of the 
freeze-out points of K meson in the z = Q plane for the central and peripheral collisions as 
in Fig. [TJ The widths of the freeze-out distributions of K meson are smaller than those of 
(j) meson because of the higher freeze-out temperature of K meson. 
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FIG. 1: Distributions of the freeze-out points of cf) meson in the z = 0 plane for central and 
peripheral collisions of Au-|-Au at = 200 GeV and Pb-|-Pb at = 2.76 TeV, plotted 

with 2000 events. 
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FIG. 2: Distributions of the freeze-out points of K meson in the z = 0 plane for the central and 
peripheral collisions as in Fig. [H plotted with 2000 events. 
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FIG. 3: (Color online) Normalized distributions of time and transverse coordinate of the cj) freeze- 
out points in the z = 0 plane for the heavy-ion collisions as in Fig. [TJ The thin dashed lines are 
for different single events and the thick solid lines are for 2000 events. 


In Fig. [3], we show the normalized distributions of time and transverse coordinate of 
(j) freeze-out points in the z = 0 plane for different single events (thin dashed lines) and 
2000 events (thick solid lines), for heavy-ion collisions of Au-|-Au and Pb-|-Pb at the RHIC 
and LHC energies, respectively. One can see that the distributions for single events are 
fluctuated and these fluctuations are smoothed in the distributions for 2000 events. 


III. BBC FUNCTIONS FOR Au + Au AT RHIC AND Pb + Pb AT LHC 


Because of the fluctuations of the space-time coordinates of freeze-out points, the BBC 
functions for single events are fluctuated. However, these fluctuations can be smoothed in 
the BBC functions averaged over event-by-event calculations for many events as, 

,, ]4ES|G.,(k,k)G.,(-k,-k) + |G„(k,-k)p] 

]fcESG»(k,k)G,,(-k,-k) 

where Ne is the total event number. 
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(a) Au+Au, 200 GeV, 0-10% 



(b) Au+Au, 200 GeV, 40-50% 



(c) Pb+Pb, 2.76 TeV, 0-10% 



(d) Pb+Pb, 2.76 TeV, 40-50% 



FIG. 4; BBC functions of (fxj) averaged over event-by-event calculations for 2000 events for central 
and peripheral collisions of Au-|-Au at y/s^N = 200 GeV and Pb+Pb at y/s^N = 2.76 TeV. 


We show in Fig. |4]the BBC functions of 00 averaged over event-by-event calculations for 
central and peripheral collisions of Au+Au at ^JsNN = 200 GeV and Pb+Pb at ^Jsnn = 2.76 
TeV. In the calculations of the BBC functions in this work, the total event number Ne for 
each kind of heavy-ion collision is 2000, and the space-time rapidity of source is taken to 
be \ri\ < 1. Unlike the oscillations appearing in the BBC functions for the hydrodynamic 
sources with Gaussian distributions of initial-energy density js], the BBC functions for the 
hydrodynamic sources with FlC vary smoothly with respect to the particle momentum. This 
is because the oscillations in the single-event BBC functions are canceled each other in the 
BBC functions averaged over event-by-event calculations for many events. From Fig. IHone 
can see that the BBC functions for peripheral collisions are larger than those for central 
collisions. Also, it can be seen that the BBC functions for collisions at the RHIC energy are 
larger than those for collisions at the LHC energy. 

In Eqs. (IT)) and (IHD, the factors and are equal to 1 and e^*^*"*, respectively, for 

ki = k, k 2 = —ki = —k. The BBC function C(k, —k) dehned with Eq. ([T]) is proportional 
to the square of the Fourier transformation of the source temporal distribution with the 
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weight -k^'^-k'j k; (^k; -k^ + ^)]’ effects of the source spatial 

distribution in the numerator and denominator in Eq. ([1]) may cancel out partially. So, the 
values of the BBC function are mainly related to the source temporal distribution, particle 
mass, mass shift, and source temperature. Because we use the same freeze-out temperature 
for 0 meson in the calculations for the collisions with different centralities and energies, the 
main reason of the dependences of the 0 BBC function on the collision centrality and energy 
(see Fig. IHfor fixed k and m*) is that the temporal distribution of the source is narrower in 
peripheral collisions and becomes wider with increasing collision energy (see Fig. [1]). 


(a) Au+Au, 200 GeV, 0-10% (b) Au+Au, 200 GeV, 40-50% 




FIG. 5: BBC functions of K^K averaged over event-by-event calculations for 2000 events for 
central and peripheral collisions of Au-|-Au at s/snn = 200 GeV and Pb-|-Pb at s/snn = 2.76 
TeV. 


We show in Fig. |5]the BBC functions of K^K~ averaged over event-by-event calculations 
for central and peripheral collisions of And-An at y/sNN = 200 GeV and Pb-I-Pb at s/snn = 
2.76 TeV. As compared to the results of the BBC functions of 00, the BBC functions of 
K^K~ are much smaller. This is mainly because the higher freeze-out temperature and 
smaller mass of K meson lead to larger values of in the calculations of BBC 

functions. 
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For anisotropic sources, the anisotropic source velocity may 


the BBC function on the direction of particle momentum 


ead to the dependence of 


29 |. We show in Fig. |6]the 


dependence of the average BBC functions of (jxp, ((^(k, —k))|k|, on the cosine of the particle 
azimuthal angle ip for the collisions at the RHIC and LHC energies. Here, n\ is taken as 
1.05 GeV, corresponding approximately to the middle between the valley and peak of the 
BBC function (see Fig. 0]), and the momentum region averaged is 0-1 GeV/c. The BBG 


functions decrease with increasing cos ip more rapidly for pern 
source transverse velocities are more anisotropic in this case js. 


leral collisions because the 


29|. 
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FIG. 6; (Color online) Dependence of the average BBC functions of (p(p on the cosine of the particle 
azimuthal angle ip for the collisions at the RHIC and LHC energies. Here, m* is taken as 1.05 GeV 
and the momentum region averaged is 0-1 GeV/c. 


We show in Fig. [7]the dependence of the average BBG functions of cpcp, (G(k, — k))|k|, on 
the particle pseudorapidity for the collisions at the RHIG and LHG energies. Here, m* is 
taken as 1.05 GeV and the momentum region averaged is 0-1 GeV/c. One can see that the 
BBG functions decrease with increasing \y\. This is because the average source longitudinal 
velocity is higher than the average source transverse velocity for hydrodynamic sources with 
Bjorken longitudinal boost invariance, and the higher l ong itudinal velocity leads to larger 


average values of e at \y\ = 1 than at ||/| = 0 [sl, l29 |. 
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FIG. 7: (Color online) Dependence of the average BBC functions of (pcf) on the particle pseudora¬ 
pidity for the collisions at the RHIC and LHC energies. Here, m* is taken as 1.05 GeV and the 
momentum region averaged is 0-1 GeV/c. 



IV. BBC FOR Cu Cu COLLISIONS AT = 200 GeV 


As we know squeezed BBC function is sensitive to the temporal distribution of the particle 
freeze-out points. The narrow temporal distribution corresponding to a small hydrodynamic 
source may lead to a large BBC. We are therefore motivated to investigate next the BBC for 
heavy-ion collisions of Cu-kCu at ^/sNN = 200 GeV, which have smaller sources compared 
to those in Au-I-Au and Pb-|-Pb collisions. 
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FIG. 8: Distributions of the freeze-out points of cj) meson in the z = 0 plane for Cu-|-Cu collisions 
at -s/snn = 200 GeV and with the impact parameter b = 0-4 and 4-8 fm. The number of event is 
2000 . 
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FIG. 9: (Color online) Normalized distributions of time and transverse coordinate of the cj) freeze- 
out points in the z = 0 plane for Cu-|-Cu collisions at = 200 GeV. The solid lines are for 

b = 0-4 fm, and the dashed lines are for b = 4-8 fm. The number of event is 2000. 


In Figs. [8](a) and [8](b), we plot the space-time distributions of the freeze-out points 
of 0 meson in the z = 0 plane for Cu-|-Cu collisions at = 200 GeV with impact 

parameter b = 0-4 and 4-8 fm, respectively. The event number for both the two impact 
parameter regions is 2000. One can see that the distributions for Cu-I-Cu collisions are 
smaller compared to the distributions in Figs. [H^a) and [11(b) for the Au-|-Au collisions. In 
Figs. Ha) andHb), we show the normalized distributions of time and transverse coordinate 
of 0 freeze-out points in the z = 0 plane for the Cu-|-Cu collisions. The normalized time and 
transverse coordinate distributions are obtained by projecting the space-time distributions 
in Fig. [8] to the time and transverse coordinate axes, respectively. One can see that the 
widths of the distributions for peripheral collisions are obviously narrower than those for 
central collisions. 

We plot in Figs. EOKa) and [TOf b) the BBC functions of 00 averaged over event-by-event 
calculations for 2000 events for Cu-|-Cu collisions with impact parameter b = 0-4 and 4-8 
fm, respectively. The BBC function for peripheral collisions is larger than that for central 
collisions. The peaks of the BBC functions reach about 1.3 and 2.5 for central and peripheral 
collisions, respectively. They are much higher than the corresponding results for Au-|-Au 
collisions at the RHIC energy and Pb-|-Pb collisions at the LHC energy, as shown in Fig. 
m The reason for this is mainly that the temporal distributions for the smaller sources of 
Cu-|-Cu collisions are very narrow. 

In Figs. CIKa) and lllf b). we show the dependences of the averaged BBC functions 00, 
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(a) Cu+Cu, 200 GeV, b=0-4 fm (b) Cu+Cu, 200 GeV, b=4-8 fm 



FIG. 10: BBC functions of (/>(/> averaged over event-by-event calculations for 2000 events for Cu-|-Cu 
collisions at -s/snn = 200 GeV and with b = 0-4 and 4-8 fm. 
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FIG. 11: (Color online) Dependences of the average BBC functions of (j)4> on the cosine of the 
particle azimuthal angle and the particle pseudorapidity for Cu-|-Cu collisions at = 200 

GeV. Here, m* is taken as 1.05 GeV and the momentum region averaged is 0-1 GeV/c. 


-(a) Cu+Cu, 200GeV - 

- (b) Cu+Cu, 200 GeV - 

-b=4-8fm - 

; --- b=o-4fm ; 

r -b=4-8fm 1 

- --- b=0-4fm ; 

. 

'.' 


((^(k, —k))|k|, on the cosine of the particle azimuthal angle 'ip and the particle pseudorapidity 
for Cu+Cu collisions, respectively. Here, m* is taken as 1.05 GeV and the momentum region 
averaged is 0-1 GeV/c. One can see that the average BBC functions almost independent 
of cosV’ for central collisions, and decrease slightly with increasing cos^’ for peripheral 


collisions because of the anisotropic transverse expansion 


29| . As discussed in section III, 


the average BBC functions for Cu+Cu collisions also decrease with the increasing \y\ like 
that for the Au+Au and Pb+Pb collisions. 
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V. SUMMARY AND CONCLUSIONS 


In the hot and dense hadronic sonrces formed in high-energy heavy-ion collisions, the 
particle interactions in medium might lead to a squeezed BBC of boson-antiboson pairs. 
In this paper, we investigate the BBC functions of 00 and K^K~ for heavy-ion collisions 
of Au-|-Au at ^/snn = 200 GeV at the RHIC and Pb-|-Pb at -s^snn = 2.76 TeV at the 
LHC, using (2 -|- l)-dimensional hydrodynamics with the FIC generated by HIJING. The 
investigations indicate that the BBG functions averaged over event-by-event calculations 
for many events for the hydrodynamic sources with the FIG are smoothed as a function of 
the particle momentum. This is different from the BBG functions for the hydrodynamic 
sources with Gaussian initial-energy distributions, which exhibit oscillations with respect to 
the particle momentum js]. The BBG functions for collisions of Au-|-Au at the RHIG energy 
are larger than those for collisions of Pb-|-Pb at the LHG energy, because the width of source 
temporal distribution increasing with increasing collision energy. The average BBG function, 
(C'(k, — k))|k|, exhibits a decrease with increasing the cosine of the particle azimuthal angle 
for peripheral collisions because of the anisotropic source transverse velocity. Based on the 
calculations, we predict that the BBG of 00 may possibly be observed in peripheral collisions 
of Au-|-Au at the RHIG energy and Pb-|-Pb at the LHG energy. The BBG of 00 is large 
for the smaller sources of Gu-|-Gu collisions at yGjviv = 200 GeV. It is of great interest to 
investigate experimentally the BBG of 00 in high-energy heavy-ion collisions. 
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